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UDC  625.122:621*. 136 
CALCULATING  T3-1PERATURE  CONDITIONS  IN  EI-!3ANKMENTS  ON  BOGS 


Moscow  TRANSPORTNOTE  STROITEL ' STVO  in  Russian  Jul  76  pp  hh-k6 

[Article  by  G.  S.  Pereselenkov,  candidate  of  technical  sciences:  "Relia- 
bility of  Calculations  of  Temperatiire  Regimes  of  Embankments  on  Bogs"] 

[Text]  The  designing  of  embankments  on  bogs  (conditions  such  as  on  the 
BAM  [Baykal -Amur  Railroad]  raises  the  problem  of  forecasting  the  tenpera- 
tiire  regime  at  the  base  and  in  the  body  of  the  embankment.  This  is  espe- 
cially important  for  segments  of  insxilar  and  high-teraperatiure  frozen  ground. 

Existing  methods  of  calculating  engineering  tasks  for  such  conditions  are 
based  on  principles  of  construction  a model  of  the  base  soils  and  schema- 
tizing actual  conditions  of  the  base's  interaction  with  the  embankment 
and  the  surro\inding  medium,  also  botmdary  conditions:  the  time  of  onset 
of  the  interaction  between  the  embankment  and  the  base;  their  initial 
temperature,  the  course  of  variation  in  the  air  and  soil  temperature; 
and  so  on. 


This  approach  to  calculations  in  and  of  itself  conceals  a certain  amount  of 
danger  of  errors  in  the  results.  In  addition,  all  calculations  of  this 
type  suffer  from  a lack  of  precise  error  criteria,  because  it  is  practically 
impossible  to  verify  them.' 

In  any  attempt  to  classify  the  sources  of  possible  errors  in  these  calcula- 
tions, four  main  groups  emerge.  Errors  of  the  model,  due  to  deviations 
of  the  actual  heat  transfer  conditions  from  those  assumed  in  the  model, 
to  differences  in  the  assimed  and  the  actual  ratio  of  conductive  and  con- 
vective heat  exchange,  and  to  inaccurate  model  reproduction  of  the  actual 
properties  of  the  embankment  earth  and  the  base. 

Errors  in  initial  data,  due  to  discrepancies  between  the  parameters  of 
actually  poured  soils  and  those  assimed  in  the  design  and  calcination, 
and  to  great  variety  and  spread  of  parameter  ’/alues  of  the  soils  of  the 
bog  itself  (See  TRANSPORTNOYE  STROITEL ' STVO  No  5,  1970,  p 39).  Errors 
in  boundary  conditions,  which  must  include  deviations  from  calculated 
values  in  the  actual  values  of  temperatures,  the  type  and  conditions  of 


The  results  of  thermotechnical  calculations  are  among  the  basic  initial 
data  in  the  designing  of  transverse  profiles  of  an  earthen  bed,  and  desig- 
nation of  the  size  of  the  "safety  factor"  to  prevent  deformations  or  to 
reduce  them  to  a level  permitting  satisfactory  operation  of  the  roadway 
does  not  always  achieve  the  purpose,  especially  in  the  projecting  of  em- 
bankments according  to  Principle  II  (maintaining  the  degree  of  freezing 
in  the  base).  For  this  reason,  it  appears  advisable  to  determine  the 
reliability  of  these  results  and  then  to  undertake  to  stipulate  the  extent 
of  their  reliability  by  norm.  This  requires  concerted  efforts  to  reduce 
the  effect  of  possible  sources  of  error  in  all  the  four  main  groups 
mentioned  above,  or  to  compensate  for  them. 

Of  considerable  interest  are  calculations  of  reliability  in  determining 
the  depth  of  season  thawing  and  forecasting  the  level  of  stabilization 
in  the  upper  boundary  of  frost  (VGM) . 

As  is  well  known,  the  depth  of  seasonal  thawing  is  determined  by  examining 
conditions  of  heat  exchange  in  soils  by  solving  the  task  of  three-dimen- 
sional thermal  process.  For  practical  engineering  calciilations , the  three- 
dimensional  task  can  be  reduced  to  a two-dimensional  or  one-dimensional 
scheme  and  solved  with  varying  degrees  of  refinement , for  which  many  years 
of  research  in  this  field  have  developed  calculation  formulas.  Recently, 
algorithms  and  programs  have  been  developed  for  making  the  necessary 
calculations  by  means  of  computers. 


Many  years  of  experience  have  established  the  sphere  of  rational  application 
of  each  of  these  calctilations  methods . Further  searches  for  new  calciilation 
formulas  constitute  a separate  problem. 

For  this  reason,  to  illustrate  the  calculation  of  reliability  in  determining 
the  depth  of  season  thawing  (freezing)  we  can  make  use  of  Stefan's 

formula,  the  simplest  and  therefore  the  most  widely  used  [l]: 

^t(f)  “ '^^^t(m)'*'av.val.^t(f ) (l) 


where  is  the  coefficient  of  heat  conductivity  of  thawed  (during  thawing) 

^ of  frozen  (during  freezing)  soil; 

^av  val  average  temperature  of  the  soil  surface  (air)  during 

the  thawing  (freezing)  period; 
is  the  duration  of  the  thawing  (freezing)  period; 


'l.h. 


is  the  amount  of  latent  heat  in  the  soil,  absorbed  by  the  ice 
during  thawing  or  given  off  by  water  during  freezing. 


Q,  . = 0.8^  W I; 

l.h.  Ysk  grav 


(2) 


where  Ygj^  is  the  volume  weight  of  the  soil  skeleton; 


W is  the  gravimetric  moisture  content  of  the  soil; 
grav 

I is  the  iciness  of  the  soil. 


The  volume  weight  of  the  soil  skeleton,  its  gravimetric  moisture  content, 
iciness,  and  the  closely-associated  coefficient  of  heat  conductivity  can 
have  various  values  not  only  within  one  bog  but  also  in  adjacent  cross 
sections.  Average  soil  surface  temperatures  and  the  duration  of  the 
thawing  and  freezing  period  can  also  vary  in  different  years.  All  this 
makes  it  possible  to  view  the  depth  of  thawing  (freezing)  as  a system 
of  random  values: 


where  is  the  random  values  of  soil  characteristics; 


(3) 


A = ^ / 


2.5X^(f)Y‘V' 

' ' sk  gray 


Q- 


i.h. 


(M 


t is  the  random  value  of  average  soil  surface  (air)  temperature  over 
the  thawing  (freezing)  period: 


T = /t  7 ; 
av.val. 

D is  the  random  value  of  the  duration  of  the  thawing  (freezing) 
period: 


(5) 


D = 


’^■’^t(f ) • 


(6) 


The  spread  of  values  for  each  of  these  random  values,  as  a rule,  can  be 
assessed  from  field  date  of  engineering-geological  siirveys,  observations 
of  meteorological  stations,  and  frozen-groxmd  studies  along  the  route. 

Taking  the  distribution  of  the  probably  values  of  these  random  values  to 
be  subject  to  the  normal  law,  which  can  be  assumed,  and  proceeding  on 
Lyapunov's  theorem,  since  the  values  of  each  of  them  depend  on  a large 
nxmiber  of  influencing  factors,  it  is  possible  to  utilize  Chebyshev's 
inequality  (the  "three  sigma"  rule)  and  find  the  standard  deviation  for 
each. 

Accordingly,  it  is  possible  to  determine  the  standard  deviations  of  the 
thawing  (freezing)  depths  as  functions  of  the  random  value  of  arguments 
making  up  the  system  [2]. 


This  distribution  of  this  function  should  also  be  subject  to  the  normal 
law  and  can  be  written  as: 


Hq)  = 


2a^H. 


^^tj  - (%.  - 


20 


202 
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(7) 


where  E(.;  represent  the  depths  of  thawing  (freezing)  as  functions  of 
random  arguments  of  soil  parameters  and  thawing  (freezing) 
periods,  respectively; 

mHy^;  mllf.;  mH^  represent  mathematical  expectations  of  corresponding 
functions  of  random  arguments; 

Ogy^;  OjjQ  represent  standard  deviations  of  corresponding  functions. 


If  we  now  take  10  percent  to  the  maximum  deviations  between  actual  thawing 
(freezing)  depths  and  calcxalated  ones,  the  reliability  of  calculation  of 
the  thawing  (freezing)  depth  will  be  expressed  by  the  probability 
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®'^^VTCal ^ eict< 


cal 


Ht-,  H^)dHAdHTdH^  = — 


I . I I . I I . I 

= J f \ /(Wa; 

0 . 9 o'!  9 o'!  9 


(8) 


where  $ is  the  normalized  Laplace  function; 

H^v  -yal  average  C't.lCTilated  depth  of  thawing,  computed  for 

the  average  values  of  A,  T,  and  D. 


For  example,  for  sandy  loam  soil  with  the  following  initial  data: 

tav  = (11°— 10°)C;  = (200—190)  days;  A^  = (l.lU—l.lo)  kcal/m-hr-deg; 

Xgjj  = (1.6— 1.2)  g/cm^;  Vgrav  * (!*0— 30)  1 = (0.8— 0.6);  = 2.01  m. 


H = 2.U3  m,  H . =1.72  (Fig.  2),  the  reliability  of  determining  the 

max  min  10 

thawing  depth  came  to  0.3. 

Increasing  the  degree  of  reliability  of  calculation  requires,  above  all, 
more  accurate  meteorological  data.  In  particular,  if  the  period  of  posi- 
tive temperatures  in  this  calculation  were  specified  within  197 — 193  days, 
the  reliability  of  calculation  would  be  increased  to  0.7. 


BCMnpu  Hcp.p  BrMnpuHmax 

(d)  (e) 


Fig.  2 
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d.  VOl  at  H, 


av . cal . 


e,  VGM  at  H 


max 


Reliability  in  predicting  the  level  of  stabilization  of  the  upper  boundary 
of  frozen  group  can  be  determined  by  representing  the  embankment  and  the 
thickness  of  the  active  layer  iinder  it  as  a technical  system  whose  mode 
of  operation  — thawing  or  freezing — varies  randomly.  "Breakdown"  of  the 
system  can  begin  in  each  mode  if  its  duration  is  greater  than  is  necessary 
for  full  thawing  or  freezing. 

The  states  of  the  system  (Fig  3)  will  be: 

1)  warming  mode;  the  system  is  "operating"  as  thawing  proceeds; 

2)  warming  mode;  the  system  "breaks  down"  as  VGM  declines; 

3)  cooling  mode;  the  system  is  "operating"  as  freezing  proceeds; 

U)  cooling  mode;  the  system  "breaks  down"  as  cold  accumulates. 


i 


Fig.  3 


The  intensity  of  the  flow  of  events  xinder  which  the  system  converts  from 
the  thawing  to  the  freezing  mode  is  proportional  to  the  ratio  of  the  heat 
released  into  the  atmosphere  to  the  amount  of  heat  entering  the  soil.  The 
reverse  transition  takes  place  under  the  effect  of  the  flow  of  events  y^, 
proportional  to  the  inverse  ratio  "Breakdown"  of  the  system  in  the  thawing 
mode  takes  place  under  the  effect  of  the  "flow  of  breakdowns"  Pg,  propor- 
tional to  the  ratio  of  the  difference  of  the  maximum  and  the  average  cal- 
culated amount  of  heat  entering  the  soil  to  the  average  calculated  amount 
of  heat  entering  the  soil.  "Breakdown"  of  the  system  in  the  freezing  mode 
y^  takes  place  with  respect  to  the  analogous  ratio  of  the  ajnount  of  heat 
escaping  into  the  atmosphere: 


Ul  = VT(-) 


yo  = t 


av. 


max  max 


(+)  - tavT(+) 


(9) 


^ '1  ~ '^avcax'^max^”)  ~ tavTC-) 


tav"^-)’ 


tavT ( - ) 


For  the  probability  of  the  operating  states  of  the  system  we  can  write 
Kolmogorov's  equation  [3]: 


dUn 


— - -{y^  + yg)?!  + y3P2 


(10) 


P-  ’ -'“3  *■  * "in- 

where  p^  and  pg  are  probabilities  of  states  1 and  3 respectively. 


J 


The  initial  conditions  in  which  these  equations  are  integrated  are  deter- 
mined by  the  initial  operatir.g  mode  of  the  system. 


For  example,  if  the  embankment  has  been  poured  in  the  summer,  the  system 
begins  to  operate  in  the  thawing  mode  with  these  initial  conditions: 


t = 0;  P2  = 1;  P2  = 0. 

Filling  out  the  initial  data  of  the  preceding  example  with  average  air 
temperatures  and  the  duration  of  the  frost  period: 

tav.val.  = -(l8°n6°)C;  t^  = (175^165)  days, 

which  makes  it  possible  to  assume: 

Ml  = l.i+;  V2  = ^3  = 

and,  as  a result  of  calculations,  to  transform  the  equations  of  (lO)  into 
the  form: 


jT  = + 0-T1P2; 

dp2 

— — = -0.8p2  + l.Up^, 


integrating  for  the  initial  mode 
t=o;  P2=l;  P2=0i 


we  get 


Pl(t)=0.337  + 0.663 

Pj(t)=0.663(e-»-°'>9U,-2'15*). 


Accordingly,  for  the  initial  mode 
t=0;  P2=0;  P2=l; 

(i.e.,  for  winter  pouring  of  the  embankment),  we  get 


Pl(t)  = 0.337(.-''-°5‘^-"-l’‘)i 

P2(t)  = 0.663Ce"'^-09t+o.337e"^'^9t. 


The  reliability  of  the  system  in  both  cases  will  be: 


p(t)  = p^(t)+p2(t)=e“0'09t. 


As  can  be  seen,  this  reliability  is  not  very  high — in  just  five  or  six 
years  the  maintenance  and  breakdown  of  the  VGM  level  will  become  equally 
probable . 

However,  an  examination  of  the  reliabilities  of  the  "operating  state"  of 
the  system  in  the  thawing  mode  pl(t)  and  in  the  freezing  mode  P2(t)  shows 
that  during  freezing  the  system  "operates"  approximately  twice  as  reliably. 
This  makes  it  possible  to  assume  that  in  this  case  there  will  be  no 
decline  in  V®I. 

At  the  same  time,  putting  the  system  "into  operation"  during  the  winter 
somewhat  increases  its  reliability  in  the  freezing  mode  (by  — 

that  is,  the  probability  of  degradation  of  the  frozen  ground  declines. 
F\irther  refinement  in  reliability  calculations  will  make  it  possible  to 
assess  the  design  of  embankments  in  bogs  in  more  detail  and  to  pred'’‘ct 
their  interaction  with  the  base  in  time. 
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